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ABSTRACT 

We calculate the spectra expected from unsteady relativistic wind models of 
gamma-ray bursts, suitable for events of arbitrary duration. The spectral energy 
distribution of the burst is calculated over photon energies spanning from eV to TeV, 



> 

OS, 

^f) , for a range of event durations and variability timescales. The relative strength of the 

emission at different wavelengths can provide valuable information on the particle 
O ' acceleration, radiation mechanisms and the possible types of models. 

Subject headings: gamma-rays: bursts - radiation mechanisms: non-thermal - shock 
^ waves - relativity 
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1. Introduction 

Relativistic wind models of Gamma Ray Bursts (GRB) are a recent development of the 
dissipative relativistic fireball scenario which is a natural consequence of observationally set 



requirements irrespective of the GRB distance scale (e.g. Meszaros 1995). The observed bursts are 
expected to be produced in optically thin shocks in the later stages of the fireball expansion. Two 
types of shocks have been considered, based on different interpretations of the burst duration and 



variability. First, the "external" deceleration shocks ( [Rees Meszaros 1992 ; Meszaros & Rees 



1993 ) that develop at r^ec, due to the unavoidable interaction of the relativistic ejecta with the 
surrounding medium, give bursts whose duration is determined (in the "impulsive" regime) by the 
relativistic dynamic timescale td ec ~ r dec/cT 2 , where V is the Lorentz factor of the expansion (e.g. 



Meszaros fc Rees 1993| ; [Katz 1994 ; Sari, Narayan & Piran 1996). Second, the "internal" dissipation 



shocks in an unsteady relativistic wind outflow ( Rees fc Meszaros 1994] [RM94]; Paczyhski & 
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Xu 1994), where both the burst duration t w and time variability t var < t w , are characteristic 
of the progenitor mechanism (e.g. a disrupted disk accretion timescale, dynamic or turbulent 
timescales, etc.). Simple "external" shocks tend to produce fairly smooth light curves with a fast 
rise followed by an exponential decay (FRED), as observed in some bursts. However, many bursts 
have multi-peaked, irregular light curves; those may find a more natural explanation (Fenimore, 
Madras & Nayakshin 1996) in the relativistic wind scenario mentioned above. 

While the spectral properties of impulsive "external" shocks have been studied in some detail 
flMeszaros, Rees fc Papathanassiou 1994| [MRP94] ; | ari et al. 1996| ), so far, only rough estimates 



have been made for wind "internal" shock spectra ( pVleszaros fc Rees 1994 ; Thompson 1994 ). It is 



of great interest to explore the spectral properties of both types of shock scenarios, particularly 



in view of the increasingly sophisticated analyses of both GRB light curves (e.g. Fenimore et al. 



1996; Norris et al. 1996) and spectral characteristics (e.g. Band et al. 1993[ ). In this letter, we 



investigate the spectral properties of "internal" shocks for representative parameter values and 
discuss those properties in the light of anticipated results from HETE. 



2. Physical Model and Shock Parameters 

In the unsteady relativistic wind model, the GRB event is caused by the release of an amount of 
energy E = 10 51 ^5i erg, inside a volume of typical dimension ri = 3 x 10 10 t mr cm, over a timescale 
t w (RM94). The average dimensionless entropy per particle r\ = (E/(Mc 2 )} — 1 = (L/(Mc 2 )} 
determines the character of the overall flow. The bulk Lorentz factor of the flow increases as 
r oc r at first and later saturates to a value r = T s at r <5 r s ~ r^r s , where T s ~ r\ for the 
values of r\ that we consider here (eq. Q). It is assumed that rj varies substantially, Arj ~ rj, on 
a timescale t var . As an example, we consider shells of matter with two different rj values ejected 
at time intervals t var . After saturation these will coast with bulk Lorentz factors ~ Qf^s and 
T r ~ q r F s , where q r < qj ~ 1. Shells of different rj catch up with each other at a dissipation 
radius « TfT r ct var 3 x 10 14 q 2 T 2 2 t var cm (see RM94, eq. [4]), where q = ^/qjql- and 
the subscript 2 (3) stands for quantities measured in units of 10 2 (10 3 ). This gives rise to a 
forward and a reverse shock which are mirror images of each other in the contact discontinuity 
(center of momentum, CoM) frame. The CoM moves in the lab frame with Tf, ~ y/T fT r , for 
T s S> 1. In the CoM frame the two shocks move in opposite directions with small Lorentz factors 
T' sh - 1 ss (1/2) [(T//r r ) + (Tr/Tf)] 1 / 2 - 1 > 1(T 2 , for T f ,T r > 1, while in the lab frame both 
shocks move forward with Lorentz factors Tf and T r respectively. The shocks compress the gas 
to a comoving frame ( primed quantities refer to the comoving frame, which is the CoM frame ) 
density 

where 9 2 is the normalized opening solid angle of the wind (1 if spherical). The bulk 
kinetic energy carried by the flow is randomized in the shocks with a mechanical efficiency 
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£sh ~ (<?/ + Qr — 2g)/(<?/ + q r ) (eq. [1] in RM94). In order for the shocks to produce a non-thermal 
spectrum and radiate a substantial portion of the wind energy, the dissipation must occur beyond 
both the wind photosphere (see section 2.2 in RM94) and the saturation radius. This restricts rj 
to the range 

33 G#H 1/5s » S8o (pnH 1/4 - (2) 

The collisionless "internal" shocks accelerate particles to relativistic energies with a 
nonthermal distribution. We parameterize the post-shock electron energy by a factor k, so that 
the average random electron Lorentz factor is 7 e ~ nV s h ~ ^ < m p /m e (with the protons' being 
7 p ~ r s /i). An injection fraction £ of the post-shock electrons is accelerated to a power-law 
distribution (dN{^) ~ 7~ p <i7, for 7 m j n < 7). We take p > 3, so that most of the energy is carried 
by the low energy part of the electron spectrum ( r y m i n ~ ((p — 2)/(p — 1))k). A measure of the 
efficiency of the transfer of energy between protons and electrons behind the shock is given by 



e pe ~ Qk,/ '(m p /m e ) ~ (K3. Here, we will use C = 1 an d ft < 10 3 ( but see Bykov & Meszaros 1996) 



and p = 3 which is consistent with the fits to observed spectra (Hanlon et al. 1994). The electrons 
get accelerated at the expense of the randomized proton energy behind the shocks, which, for an 
individual shock (of which there are t w /t var ), represents a fraction (t var /t w )e pe £ s h of the total 
energy Eq. 

Magnetic fields in the shocks can be due to a frozen-in component from the progenitor, or 
may build up by turbulence behind the shocks. We parameterize the strength of the magnetic 
field in the shocks by the fraction A of the magnetic energy density to the particle random energy 
density (u B = \rirripC 2 ). The comoving magnetic field therefore is 



1.72 x 10 4 / E 51 e sh A 

~ q *ri V en w ti ar m v S2 ^ {6 > 

The relativistic electrons will lose energy due to synchrotron radiation, and inverse 
Compton (IC) scattering of this radiation. The respective radiative efficiencies are 
£ sy = t'sy 1 /^'- 1 + t'r 1 + t'- 1 ), and e ic = t'r 1 /^ 1 + tr} + t'' 1 ) . The timescales are 
defined below. The comoving synchrotron timescale is determined by the least energetic electrons : 



a2 f t 2 noT 5 

L w b var I* s 



t' sy « 8 x 10 8 / lmm B' 2 = 5.2 ^^Hl 'J^Jl ms , (4 ) 



where the subscript sy (ic) refers to the synchrotron (IC) emitting shock. The IC cooling depends 
on, and competes with, synchrotron cooling. IC cooling dominates if the magnetic field is weak, 
a large fraction of electrons are accelerated and share the protons' momentum very efficiently 
(i.e. A/(£Kj C3 ) <C 1); while synchrotron cooling dominates in the opposite case. The IC timescale 
(t' ic «3x 10 7 /u'syjmin) for the two limiting cases is 
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4 nl, +2 „. r .-. . ...... 

ms if IC dominates 



(e sh E 51 < 7| > «, 



»'3 



(5) 

4.95 — S2 ms if synchr. dominates. 

£ sh E 51 J CXk . c3 < 7 2 > 



Here 7* is the electron Lorenz factor that corresponds to the peak emission frequency, 

7* = max [7mm, Tofts]) labs is defined below equation (|) and < 7 | >= 2 x 10 _6 7^ n In ( r y m ax/l*)- 

The lab frame shell width (shocked region) is Ar ~ r d^ b 2 ~ ctct var , and the comoving 
crossing time t' ex « 10 2 (ag 2 ) r s2 t var s provides an estimate for the adiabatic loss time. 

The spectrum of an "internal" shock burst consists of two synchrotron and four IC components 
(coming from all shocks' combinations). The synchrotron components are characterized by up to 
three break frequencies, given below in the lab frame: 

i) The 7 m i n break frequency z/ = 8 x 10 5 B ,r y~ nin , which in the lab frame gets blue-shifted by 
Tf (T r ) for the forward (reverse) shock. Using = [qf,q r ] to refer to the shock, we have 



ii) The self absorption frequency v abs . If the electron power law starts at low enough energies, 
the radiation field becomes optically thick at a frequency determined by \m e ^ abs v^ hs = F' v b , and 
is obtained by solving a non-linear algebraic equation; in the limit v abs S> v m in it is 

hu abs « 18.4 -L ' (C^ 3 ) 2/7 eV for 7^ < Tofts, (7) 

where -y abs = 9 x 10 3 yj Uab s / {q sh T S2 B) . 

iii) The frequency where the photon spectrum of the minimum energy electrons becomes 
optically thick (v RJ ), and below which it assumes a Rayleigh- Jeans spectrum slope. This happens 
at u min when j min < -y abs ; if j min > j abs it is 

, - 1t . qsh ( E 51 e sh \^ f V2 \ 2 / 5 fCe sh \ 3/5 1 v 



iv) A cutoff is expected at v max = {^ ma x hmin) 2 v m in, where j max is the electron energy in each 
shock at which the electron power law cuts off due to radiative or adiabatic losses. It is determined 
by tacci'fmax) ~ t' coo i (imax ) , where t acc is the electron acceleration timescale in the shocks (a 
multiple 10 x Aiq of the inverse gyro-synchrotron frequency), t acc ~ 3.57 x 10~ 6 AiqB' (r^)^ 1 ^/ 
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s, and t coo i is the minimum of all the radiation and adiabatic cooling timescales involved, 

Each synchrotron component can be characterized by three frequencies in ascending order, 
v % S y j , where i = 1, .., 3 and j = 1 for the reverse and j = 2 for the forward shock. Similarly, the pure 
and combined IC spectra are characterized by the frequencies v\ c ^ « (4/3) K% c v l sy ; where, j = 1, ..,4 
(1 corresponds to IC reverse ,2 to IC forward, 3 to IC reverse-forward and 4 to IC forward-reverse), 
and if j is odd, I = 1, otherwise I = 2. The shape of each component depends on the relationship 
between the relevant ^mimlabs and Imax- In a power per logarithmic frequency interval plot, the 
fluence of each component exhibits a peak of Si = 1.6 x 1CP 6 (£51 / (9 Z?2s) 2 ) £sh£i(Ki 3 erg/cm 2 , 
where i = 1(2) for synchrotron (IC), and -D28 is the luminosity distance corresponding to z ~ 1 
for a flat Universe, with i7 Q w 80 (Dx = (2c/H )(l + z — \f\ + i 3 Gpc)). The spectrum is 
obtained by adding up these six components. In practice, the forward and reverse components 
have values very close to each other and they essentially merge. The resultant spectrum is then 
checked for the effects of pair production; the 77 optical depth is calculated for each comoving 
frequency above m e c 2 /k using the number density of photons above the corresponding threshold, 
as obtained from the initial spectrum; finally the spectrum is modified accordingly. We note that 
most of the scattering in our spectra occurs in the Thomson regime, unlike in [Sari et al. 1996 , 
who consider large 7 in the framework of impulsive shock models. ( Klein-Nishina corrections may 
become relevant in a few k 10 2 cases, but at frequencies which lie in the 77 absorbed part of the 
spectrum) . 

3. Typical Wind Spectra 

We discuss here the properties of some representative spectra. We assume a total event 
energy of E = 10 51 erg and a geometry of a spherical section, or jet, of opening angle 8 = 0.1 (the 
physics is the same as in a spherical wind , provided 6 > T" 1 ). We have investigated a range of 
dynamic parameters (77, t w and t var ), and used qj = 2,q r = 0.5, and kj = n r . 

In figure [l] we present spectra for a long burst (t w = 100 s) and in figure |2| for a short 
one (t w = 1 s), for different 77, A and k values. The sharper spectral features present would be 
smoothed by inhomogeneities in a real flow. 

The range of 77 where a substantial fraction of the wind energy is radiated with nonthermal 
spectra is considerably restricted (eq. Q) and implies lower values than in impulsive models. 
Larger 77 may be difficult to produce in a natural way, while lower values lead to shocks below 
the photosphere that make bursts too dim to be observed. The allowed volume of parameter 
space is not very large ( [Papathanassiou Meszaros 1996| [PM96]). Nonetheless, it allows for an 
appreciable variety of spectra, since models differing only slightly in 77 can produce significantly 
different spectra. This is due to the strong dependence of the photon number density on T s 
(i.e. n' 7 oc r~ 7 ) which determines pair production. Lower k spectra are less affected by pair 
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production, since they contain less energetic electrons and therefore fewer photons. Generally, 
other parameters being equal, the higher r\ models tend to produce spectra that span over a wider 
range of frequencies. For spectra like the majority of the ones observed by BATSE ( Fishman & 



Meegan 1995| ) values of k < 10 2 and A > 10 1 are excluded. Low k spectra are either too dim 



or too soft; high A brings up the synchrotron component and may be appropriate for a small 
percentage of bursts with low frequency excess ( Preece et al. 1996| ). In general, for a given k value 



a wide range of A values (3-6 orders of magnitude) is allowed, the trend being that higher k 
values must combine with lower A values, fairly independently of r] (within the allowed range of 
values). This trend is due to the fact that lower magnetic fields require higher electron energies in 
order for the break to fall in the BATSE window (e.g. eq. 0, ). High energy power laws (like 
those reported in Hanlon et al. 19"94| ) are common (see fig. || and right column of fig. |]). For a 



more complete discussion see PM96. 

The effect of a longer t w is to push the pair cutoff to higher frequencies, because, for 
fixed Eq, the flux and photon density are lower. For the same rj, k and t var , longer bursts 
require higher A. A greater t var with the rest of the parameters unchanged would again require 
higher A in order to produce observed-like spectra. The cases considered here were chosen with 
tw < t d ec ~ (E /n ex m p c 5 r] 8 ) 1 / 3 s. 



4. Conclusions 

An unsteady relativistic wind provides an attractive scenario for the generation of GRB, since 
it requires smaller bulk Lorentz factors than the impulsive models, i.e. it can accommodate higher 
baryon loads (RM94). In addition, the lack of kinematic restrictions (e.g. Fenimore et al. 1996| ) 



allows it, in principle, to describe events with arbitrarily complex light curves. We have calculated 
spectra from optical through TeV frequencies for bursts with a range of durations and variability 
timescales. At cosmological distances the total energies and photon densities implied by the model 
are likely to turn those spectra optically thick to 77 pair production. Most of them may therefore 
be missed by, or show a high energy cut-off in, the EGRET window, but they would be prominent 
in the BATSE and HETE gamma-ray windows. Low frequency (down to 5 keV) excess reported 
recently ( |Preece et al. 1996 ) may be attributed to a pronounced synchrotron component due to a 



relatively high magnetic field, in a fraction of the bursts. 

The continued propagation of an unsteady wind flow should eventually lead to its deceleration 
by the surrounding medium. If the latter is of any appreciable density, it could lead to another 



burst with the "external" shock characteristics (MRP94, Meszaros fc Rees 1994 ), provided that 
not all of the wind energy was radiated away by the "internal" ones (e s /j < 1). The spectra of the 
"internal" shocks are different from those coming from the "external" ones, the main differences 
being that the former cover a narrower range of frequencies and most have high energy cutoffs 
due to pair production opacity (PM96). If GRB progenitors have escaped their parent galaxies 
and are in a low density intergalactic medium (n ex < 10~ 3 cm -3 ), the "external" impulsive shocks 
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would be of long duration (> 3 x lO 3 //^ 3 s ) an d low intensity, and most would be totally missed, 
hence the GRB could be entirely due to unsteady wind "internal" shocks. If the GRB occur in a 
denser medium (e.g. the galactic ISM), the last peak of multiple-peaked GRB would be smooth 
and FRED-like but the earlier peaks, due to unsteady wind "internal" shocks, may be arbitrarily 



complicated. Those might also be responsible for delayed GeV emission ( Meszaros Rees 1994 ). 



A number of unknown parameters enter into the calculation of GRB models, and HETE may 
provide the information required to narrow the range of their allowed values, as well as a test 
for the general scenario of unsteady relativistic winds. For the HETE sensitivities indicated in 
figures |l| and ^ we expect that some bursts will show a simultaneous X-ray counterpart. However, 
detection by the Ultraviolet Transient Camera (UTC) on HETE should be rare for the majority 
of bursts (note though that figures [l] and || refer to z ~ 1 distances, and that at a few x 10 2 Mpc a 
bright burst would be 30 — 100 times brighter, thus increasing the likelihood of detection by the 
UTC). For the majority of faint bursts, a detection in the ultraviolet would be possible only for 
low values of k, resulting in a very broad and fiat spectrum with upper cutoffs, if any, only at the 
highest energies. 
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Fig. 1. — Spectra for E§\ = 1,6 = 0.1, t w = 100s, t var = 40ms, 77 = 63 (left column), r\ w 160 ( right 
column), and for different values of the magnetic field and the electron acceleration efficiency. The 
long dashed line shows the approximate threshold and window for GRO's BATSE and EGRET 
experiments, while the dotted-dashed one is for all the experiments on board HETE. 
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Fig. 2. — Spectra for = 1,8 = 0.1, t w = ls,t var = 40ms, rj = 200 and for representative 
values of the magnetic field (from bottom to top log{\) = 0, —4, —8, —12) and electron acceleration 
efficiency. GRO and HETE windows are included. 



